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ABSTRACT 

A path  finding  method  and  a stochastic  time  integration  scheme  for  the  simulation  of  thermally 
activated  magnetization  processes  are  introduced.  The  minimum  energy  path  and  the  saddle 
points  for  the  thermally  induced  transitions  between  the  ground  states  of  NiFe  magnetic  nano- 
elements are  calculated. 

INTRODUCTION 

With  decreasing  size  of  magnetic  nanostructures  thermal  effects  become  increasingly 
important.  Prominent  examples  are  magnetization  noise  in  magnetic  sensor  elements  [1-3]  and 
the  thermal  stability  of  magnetic  MRAM  (Magnetic  Random  Access  Memory)  cells  [4]  or 
magnetic  storage  media  [5].  Magnetic  sensors  require  a high  sensitivity  so  that  small  magnetic 
fields  can  be  detected.  On  the  other  hand  thermal  fluctuations  which  will  lead  to  thermal  noise 
should  be  suppressed.  The  free  layer  of  a multilayer  sensor  element  is  soft  magnetic  and  may 
have  a size  well  below  one  micrometer.  Thermally  induced  magnetization  processes  may  cause 
local  or  global  magnetization  rotations  which  cause  the  magnetization  noise.  With  decreasing 
lateral  extension  of  the  elements  thermal  fluctuations  become  more  pronounced.  Magnetic 
storage  elements  require  a low  and  well  defined  switching  field  which  in  practice  is  limited  by 
the  current  through  the  write  line  in  an  array  of  MRAM  cells.  On  the  other  hand  the  shape  or  the 
induced  anisotropy  should  guarantee  a life  time  of  a stored  bit  of  about  10  years.  Again  the 
energy  barrier  for  thermally  induced  magnetization  reversal  decreases  with  increasing  size  of  the 
storage  elements.  The  corresponding  time  scales  differ  by  several  orders  of  magnitude:  Thermal 
noise  arises  on  a time  scale  of  a few  nanoseconds;  thermally  induced  switching  of  the 
magnetization  over  energy  barriers  extends  over  seconds  to  years.  Random  thermal  fluctuations 
of  the  magnetization  are  the  underlying  physical  process  which  causes  both  thermal  noise  and 
spontaneous  switching.  The  stochastic  fluctuations  arise  from  the  interplay  between  the  lattice 
vibrations  and  the  magnetization. 

A micromagnetic  system  will  be  close  to  a local  minimum  of  the  total  magnetic  Gibbs’  free 
energy.  Thermal  fluctuations  of  the  magnetization  cause  the  magnetization  to  wander  around  this 
minimum.  Occasionally  the  system  will  reach  a region  next  to  a saddle  point.  The  system  may 
cross  the  energy  barrier  and  move  into  the  basin  of  attraction  of  a different  energy  minimum. 
This  process  can  be  described  by  the  Neel-Brown  theory  [6,7].  The  relaxation  time, 

T = /„■'  exp(£b  I kBT) , is  the  inverse  of  probability  per  unit  time  for  crossing  the  barrier  Eb.  The 
attempt  frequency,/0,  depends  on  material  parameters,  like  anisotropy,  particle  shape,  and 
damping  [8],  Its  value,  which  ranges  ffom/o  = 109  Hz  to  fo  = 10i2  Hz,  sets  the  time  scale  for 
thermally  assisted  magnetization  reversal,  r0  = = 1 ns . 

The  theoretical  treatment  of  thermally  induced  magnetization  processes  starts  from  the 
stochastic  Landau-Lifshitz  Gilbert  equation  and  the  corresponding  Fokker-Planck  equation  [7], 
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The  energy  barrier  can  be  calculated  for  coherent  rotation  in  single  domain  particles  and  the 
formation  of  reversed  domains  in  thin  ferromagnetic  wires.  The  attempt  frequency  can  be 
estimated  by  solving  the  Fokker-Planck  equation  numerically  or  analytically  [7-10]. 
Alternatively,  the  stochastic  Landau-Lifshitz  Gilbert  equation  can  be  solved  numerically  for 
short  time  scales  and  small  systems  [ 1 1-14],  Recently,  numerical  solutions  of  the  stochastic 
Landau-Lifshitz  Gilbert  equation  were  reported  for  extended  micromagnetic  systems  [3,14],  Zhu 
[3]  analyzed  the  magnetization  noise  in  submicron  sized  sensor  elements  based  on  the  numerical 
solution  of  the  stochastic  Landau-Lifshitz  Gilbert  equation.  The  time  integration  of  the  stochastic 
Landau-Lifshitz  equation  is  restricted  to  small  time  scales  and  thus  is  a proper  tool  to  analyze 
magnetization  noise.  The  calculation  of  long  term  thermal  effects  needs  a detailed 
characterization  of  the  energy  landscape  along  the  most  probable  path  which  is  taken  by  the 
system  to  cross  the  energy  barrier.  Berkov  [15]  calculated  the  transition  path  of  interacting  single 
domain  particles  by  minimizing  the  action  along  the  path.  He  showed  that  a direct  minimization 
of  the  action  may  also  give  paths  through  local  maxima  which  have  to  be  excluded.  Ren  [ 1 6] 
proposed  an  elastic  band  method  to  calculate  a minimum  energy  path  in  micromagnetic  systems. 
Starting  from  an  initial  guess  for  the  path  which  connects  two  local  minima  of  the  system,  a 
highly  probable  path  is  found  by  moving  the  points  along  the  path  according  to  an  algorithm 
which  resembles  tensioning  an  elastic  band  across  a mountain.  Variants  of  elastic  band  methods 
are  commonly  used  to  calculate  transition  rates  in  physical  chemistry  [17],  Dittrich  an  co- 
workers [19]  originally  applied  the  elastic  band  method  to  calculate  energy  barriers  and  saddle 
points  in  complex  micromagnetic  systems  like  discrete  perpendicular  recording  media. 

In  this  work  we  apply  both  stochastic  time  integration  and  path  finding  techniques  in  the 
framework  of  the  finite  clement  method,  in  order  to  simulate  thermal  effects  in  magnetic 
nanostructures.  Thus  it  is  possible  to  take  into  account  complex  geometries  and  realistic  element 
shapes.  Both  methods  are  complementary.  The  stochastic  time  integration  is  restricted  to 
simulation  times  of  about  10  ns.  As  a consequence  the  calculation  of  barrier  crossing  by 

j stochastic  time  integration  is  limited  to  small  energy  barriers.  The  transition  rate  for  large 

barriers  can  be  estimated  from  the  barrier  height  which  can  be  calculated  from  the  minimum 
energy  path.  In  addition  to  the  energy  barrier,  the  elastic  band  method  provides  a global  view  of 
the  energy  landscape  such  as  local  minima  and  saddle  points  along  the  path.  The  magnetization 
processes  as  computed  from  the  stochastic  time  integration  method  and  the  minimum  energy 
path  are  compared  for  transitions  between  different  ground  states  in  magnetic  nano-elements. 

j MICROMAGNETIC  AND  NUMERICAL  BACKGROUND 

j The  micromagnetic  description  of  the  system  starts  from  the  total  magnetic  Gibbs'  free 

energy  [19] 

£ = + AT, [l  - (u  • a)J-  y (u  • H J-  7 , (u  • H„,  )jdV . (1) 

E is  the  sum  of  the  exchange  energy,  the  anisotropy  energy,  the  stray  field  energy,  and  the 
Zeeman  energy;  u denotes  the  unit  vector  parallel  to  the  magnetization,  A is  the  exchange 
constant,  Kt  is  the  uniaxial  magnetocrystalline  anisotropy  constant,  and  /» is  the  spontaneous 
magnetic  polarization.  The  integral  (1)  is  over  the  total  volume  of  the  magnetic  particles.  In  a 
stationary  state  the  magnetic  system  occupies  a local  minimum  of  ( l ).  Owing  to  thermal 
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activation  the  system  may  overcome  an  energy  barrier  and  spontaneously  move  towards  a 
different  local  minimum  of  the  energy. 

We  use  the  finite  element  method  to  evaluate  E for  complex  magnetic  systems.  The 
direction  cosines  of  the  magnetization,  «t,  are  interpolated  by  piecewise  linear  functions  on  a 
tetrahedral  finite  element  mesh.  In  order  to  calculate  the  magnetic  stray  field,  Hs,  we  use  a hybrid 
finite  element  / boundary  element  method  [20],  The  simulation  of  the  time  evolution  of  the 
magnetization  requires  to  calculate  the  effective  field,  H^,  defined  by  the  negative  variational 
derivative  of  E.  The  effective  field  on  the  nodes  of  the  finite  element  mesh  may  be  approximated 
using  a box  scheme 
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1 dE 
Mo  9m'  ’ 
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where  m!  is  the  magnetic  moment  associated  with  node  / of  the  finite  element  mesh.  The 
stochastic  Landau-Lifshitz  equation  is  a system  of  3N  Langevin  equations  with  multiplicative 
noise 
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where  the  indices  ij,k  run  over  the  three  space  directions,  and  the  index  / = 1 , . . . JV  runs  over  the 
number  of  nodes.  The  drift  term,  A I , is  the  right  hand  side  of  the  deterministic  Landau-Lifshitz- 
Gilbert  equation,  a is  the  Gilbert  damping  constant  and  y is  the  gyromagnetic  ratio.  H * is  the 
random  thermal  field.  The  thermal  field  is  assumed  to  be  a Gaussian  random  process  with  the 
following  statistical  properties: 


(tf ,U'))  = 0 and  (wl,(r),W‘./(f))  = 2 DS^Sit-f)  (6) 

The  average  of  the  thermal  field,  taken  over  different  realizations,  vanishes  in  each  direction 
i in  space.  The  thermal  field  is  uncorrclated  in  time  and  space.  The  strength  of  the  thermal 
fluctuations  follows  from  the  fluctuation-dissipation  theorem  [11]: 


D = 


akj 

|tKK 


(7) 


Numerically  the  equations  (3)  to  (5)  are  solved  using  a semi-implicit  time  integration 
method  [21].  The  time  is  divided  into  a lattice  of  discrete  points  t„.  The  Langevin  equation  is 
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solved  in  the  time  interval  t„  to  f„, , with  the  initial  condition  u‘  = «'(/„).  The  right  hand  side  of 
(3)  is  evaluated  at  the  midpoint  = (r„  +t,.,)/2  = r„  + At/2.  The  update  of  the  direction  cosines 
is  given  by 

= A' a!",...)  + (8) 

k 

with  the  noise  integral 

aw;,  = ]h'm  m.  (9) 

i„ 

The  magnetization  directions  at  the  midpoint  of  the  time  interval  is 

))/2  = k,'(0  + A«/2,  (10) 


assuming  linearization  of  uj  within  the  time  interval.  The  midpoint  value.  «/ , is  given  implicitly 
by  the  equations  (8)  and  (10)  which  is  solved  by  functional  iteration.  Numerical  tests  show  that 
about  five  iterations  are  sufficient  to  gain  the  required  accuracy  of  10  \ From  the  midpoint  value 
the  magnetization  at  the  time  r«i  is  calculated  from  (10)  and  then  normalized: 

“/(',o)=2 «/-«,'(»„).  u'=n/|Sj.  (II) 

A common  feature  of  path  finding  methods  is  the  discrete  representation  of  the  path 
connecting  the  inital  state  of  the  system  with  its  final  state.  In  micromagnetics  we  represent  the 
magnetic  states  of  a system  by  the  set  of  magnetic  moments  M = ( . . . ,m' 1 , m\  1 ,. . . | .The 

index  1=  1 runs  overall  nodes  of  the  finite  element  mesh.  For  the  configuration  space  we 
use  polar  coordinates  of  dimension  2N.  First  we  construct  a sequence  of  magnetic  states  in  such  a 
way  as  to  form  a discrete  representation  of  a path  from  the  initial  magnetization  state,  M(”,  to  the 
final  magnetization  state,  Mtn.  An  optimization  algorithm  is  then  applied  until  at  any  point  along 
the  path  the  gradient  of  the  energy  is  only  pointing  along  the  path.  This  path  is  called  minimum 
energy  path  which  means  that  the  energy  is  stationary  for  any  degree  of  freedom  perpedicular  to 
the  path.  The  minimum  energy  path  typically  represents  the  path  with  the  greatest  statistical 
weight.  From  this  path  statistical  quantities  as  for  example  transition  rates  for  the  thermally 
induced  magnetization  reversal  can  be  estimated. 

Henkelman  and  Jdnsson  proposed  the  nugded  elastic  band  method  to  calculate  minimum 
energy  paths  [17],  We  represent  a path  by  a sequence  of  images.  An  intitial  path  is  assumed 
which  connects  the  initial  magnetization  state  M10  = M11’  with  the  final  magnetization  state 
M(fl  = M'").  The  index  k runs  from  1 to  m.  The  path  is  optimal,  if  for  any  image  Mrt>  the  gradient 
of  the  energy  is  only  pointing  along  the  path  or  in  other  words  the  component  of  the  energy 
gradient  normal  to  the  path,  D,  is  zero.  If  t denotes  the  unit  tangent  vector  along  the  path,  a 
minimum  energy  path  has  the  following  property 


D,4|=V£(M'H)-[V£(M“,)  t]t=0,  for*=l m 
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The  optimal  path  can  be  found  using  an  iterative  scheme.  In  each  iteration  step  the  images 
move  towards  lower  energy  in  a direction  perpendicular  to  the  path.  So  the  image  M®  is  moved 
into  the  direction,  This  iterative  scheme  is  numerically  very  ineffective.  Thus  instead  we 
solve  a system  of  ordinary  differential  equations 

= -Dl°  fork  = 2,...,  (m-  1)  (13) 

3r 

using  an  implicit,  variable  order,  variable  time  step  time  integration  method  [22].  Here  the  time  t 
is  introduced  for  numerical  convenience  and  has  no  physical  meaning.  In  order  to  keep  an  equal 
distance  between  successive  images  a spring  force  may  be  introduced  [17]. 

RESULTS 

Small  rectangular  NiFe  elements  show  two  distinct  ground:  The  S-state  and  the  C-state. 
Thermal  fluctuations  may  induce  a transition  between  the  different  states.  The  C-state  has  a 
slightly  lower  energy  than  the  S-state.  Figure  1 shows  the  magnetic  states  during  the  thermally 
induced  switching  from  the  S-state  to  the  C-state.  The  magnetization  distributions  are  calculated 
from  the  numerical  solution  of  the  stochastic  Landau-Lifshitz  Gilbert  equation  at  a temperature 
of  T = 350K.  The  sequence  of  magnetization  configurations  in  figure  1 covers  a time  of  5 ns. 

The  transition  between  the  grounds  states  were  also  calculated  using  the  path  finding 
method.  The  initial  path  was  constructed  connecting  the  two  possible  S-states.  Along  the  initial 
path  the  different  images  are  obtained  by  linear  interpolation.  Then  a minimum  energy  path  is 
found  solving  equation  (13).  Figure  2 shows  the  energy  along  the  final  path.  In  addition  the 
magnetization  states  along  the  path  are  shown.  The  path  finding  algorithm  identifies  the  C-state 
as  a global  energy  minimum.  The  C-state  is  separated  from  the  S states  by  saddle  points.  The 
comparison  of  the  magnetization  states  presented  in  figure  1 and  in  figure  2 clearly  shows  that 
the  magnetic  states  along  the  minimum  energy  path  correspond  to  the  magnetic  states  which  are 
visited  during  the  stochastic  motion  of  the  magnetization  according  to  the  Langevin  equation. 
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Figure  1.  Magnetization  states  visited  during  the  solution  of  the  stochastic  Landau-Lifshitz 
Gilbert  equation.  The  NiFe  element  has  an  extension  of  150  x 100  x 5 nm5.  An  induced 
anisotropy  of  K\  = 500  J/m3  was  assumed  parallel  to  the  long  axis.  The  grey  scale  maps  the 
magnetization  component  parallel  to  the  short  axis  of  the  element. 
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Figure  2.  Magnetization  configurations  and  energy  along  the  minimum  energy  path.  The 
temperature  is  T=  350K.  The  saddle  point  with  the  open  symbols  is  not  shown  on  the  left  hand 
side. 
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